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ABSTRACT: We study the structure of the phase space in Hoféava-Lifshitz theory. With the 
constraints derived from the action, the phase space is described by five fields, thus there is 
a lack of canonical structure. The Poisson brackets of the Hamiltonian density do not form 
a closed structure, resulting in many new constraints. Taking these new constraints into 
account, it appears that there is no degree of freedom left, or the phase space is reduced 
to one with an odd number of fields. 
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1. Introduction 


Hořava recently proposed a gravity theory with asymmetry between time and space ia 2 
This theory is non-relativistic in the UV limit, thus it is hoped that it is UV finite. It is 
similar to a scalar field theory of Lifshitz El in which the time dimension has weight 3 if 
a space dimension has weight 1, thus this theory is called Hořava-Lifshitz gravity. Hořava 
argued that it is superficially renormalizable based on power counting and may flow to 
Einstein’s general relativity in the IR region. This work has stirred up a surge of research 
on possible applications of the theory to cosmology and black hole physics. We have no 
intention to be complete in offering the literature, for those interested in a list of these 
papers, we refer to a most recent paper pl. 

When any theory claims to be a renormalizable field theory of gravity, one must exercise 
great care, especially with a theory without general covariance to begin with, since we know 
that general covariance is hard to avoid in a theory with a massless spin 2 particle. One may 
already poses questions at face value. If a theory can flow to Einstein’s theory only when 
a cosmological constant is introduced, how can one avoid this cosmological constant?! If 
one fine-tunes parameter A in the kinetic term in the Hořava action to 1/3 thus makes the 
cosmological constant vanish, then how can flow to 1 in Einstein theory? And, how a field 
theory of gravity can explain the fact that the maximal entropy of a region is proportional 
to the area of the surface surrounding it? 

With the above questions in mind, we begin a study on constraints in this theory. 
One easily comes to doubt that whether the system of constraints of the theory makes 
sense, since not all constraints correspond to local symmetries. Hořava retains all diffeo- 
morphism symmetries in space, but gives up on time local symmetry. In a constrained 
system, constraints are normally generators of symmetries, thus they are guaranteed to 
form a closed system under the Poisson bracket. Now, the constraints corresponding to 
the lapse function have no corresponding symmetries, it is natural that they will generate 
new constraints under the Poisson bracket. Indeed they do, as we will show shortly. 

One may simply choose not to impose these new constraints. If so, then there are four 
families of constraints: one corresponding to the lapse function N, three corresponding 


1Some papers also pointed out the problem related to cosmological constant p, i. 


to shift functions N’. In the ADM canonical formalism (especially suitable for the non- 
relativistic theory of Hořava), there are 12 fields on the phase space, 6 are the spatial 
metric components gij, 6 are their canonical momenta. Upon imposing four constraints 
and three gauge symmetries, the phase space is described by five functions, and there is no 
symplectic structure on this phase space. 

Thus we need to impose new constraints derived from the Poisson brackets, but we 
shall see that there are too many new constraints, thus it appears that no degree of freedom 
is left, or the phase space is reduced to a smaller one still described by an odd number of 
fields, this seems to us a fatal problem in Hořava’s theory. 


2. Non-closure of constraints algebra and New constraints 


We start with the relativistic metric gay in the usual ADM decomposition, 


( —N?4N;N* N; ) . (2.1) 


Ni Jij 


for any theory of metric to be a theory of gravity, the lapse function N must be a function of 
both space and time, since Newtonian potential is embedded in it. Even though Newtonian 
potential can be included in the shift functions by choosing certain special gauge, to obtain 
Newtonian equation which determines Newtonian potential, the lapse function N must be 
a function of both space and time. Without Newtonian equation as a local constraint, 
many unphysical solutions will emerge. In UV region the action in Hořava-Lifshitz theory 
takes the following form 


-18 3 SK = 2) Oo 
$= f dtd?x./gN {(KijK — AK?) TCC , (2.2) 


where for simplicity we consider the first action proposed in Bl in this section, we postpone 
a discussion of the second action (which may be generated by flowing the above action 
and flow to Einstein action in the IR region) in [2] to sect.3, where similar calculations are 
carried out. The first two terms in the above action comprise the kinetic term and the 
last term is the potential term. This action satisfies the detailed balance principle, since 
CË = éy (RI — IRÖ ) is the Cotton tensor and can be obtained from the variation of a 
3-dimensional action W [g;;] P] 


A (2.3) 


where W [g;;] is the 3-dimensional Chern-Simons action. Kj; = ow (dij — ViN; — V; Ni) are 
components of the extrinsic curvature of the 3-dimensional hyper-surface with constant t 
and K = gř K;j is its trace. kw is a constant in this theory with mass dimension 1. 

In the following, we focus on the case of A = 1 so that the kinetic term in the action 
looks like that appearing in Einstein-Hilbert action. Our calculation demonstrates that 
different value of A will not change the final results qualitatively, so we can make this 
choice. Hereafter, we will adopt the units in which 167G = 1, for convenience. In the 


Hamiltonian formalism, we first compute the conjugate momenta of N, N; and gij denoted 
by 7,7’ and n, respectively. They have the explicit form 


o (2.4) 
ôN 
Pa (2.5) 
Ni 
5 6S 7 5 
r” = — = yg(K” — g” K). (2.6) 
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The basic Poisson brackets of the canonical variables are 


(NG), i)e = P- y), 
{Ni(e), 7(y)}p» = HPE- y), 
2) 


{ou (a), T"! U) = Z(E + iy), (2.9) 
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where we have adopted the convention of Poisson brackets between canonical variables 
used by 4. One should be aware of the fact that m;j is not a tensor under coordinate 
transformation, but behaves as a tensor density. In other words, 7;;/,/g is a tensor. Then 
after performing a Legendré transformation, the Hamiltonian can be derived as 


H= | wiggd’x—L= | Px(NH + NH), (2.10) 


with H and H; given by 


l 1 a 

H = —rG, in" + —./9C;;0%, 2.11 
gr Guim tga v20 (2.11) 

Hi = —2g0;09 — (20k 943 — Oiggn)T™. (2.12) 


where Gijki = 5 (Gik Gil + Gil9jk — 9ijGk1) is the inverse of de Witt metric. Since Eqs.(2.4) 
and (2.5) show that the canonical momenta of N and N’ always vanish, the two equations 
below are always satisfied 


a = {r, H} = H, (2.13) 
i = {rt, Hyp = Hi, (2.14) 
provided the following constraints 
H = 0, (2.15) 
Hi = 0. (2.16) 


Above equations are usually called as the Hamiltonian constraint and momentum con- 
straint, respectively. Once the Poisson bracket is added into the system as a part of 
structure, consistency requires that the Poisson brackets among the constraints only gen- 
erate constraints. As discussed in the introduction, we anticipate that new constraints 


are generated by the Poisson brackets among Hamiltonian constraints. To compute the 
Poisson brackets of each pair of constraints, it is helpful to first compute 


(f dx’ CHp, gijp = —C*Ongig — Gjr0i6" — Gind;C*, (2.17) 
(f Bx'C Hy, eS pp = —Og (a4 C*) + ORC + ORC, (2.18) 
which reveal the role of H;’s as generators of 3-dimensional coordinate transformation, and 


(| d°x'nH, gij}pp = = Opes. — gijgk yr"! = —2nKi;. (2.19) 


vI 


Then one obtains the following two Poisson brackets straightforwardly, 
(| axcini, | Eyan = | xia - GOA (2.20) 
{ i ax Ts J dynH}pp = f PxC ANH, (2.21) 


Eq. (2.20) tells us that H;’s form a Lie algebra corresponding to 3-dimensional diffeomor- 
phism group. Eq. (@-21) is a reflection of the fact that H behaves as a scalar density under 
this transformation. Actually, one can expect these results, since H;’s take the same form 
as in general relativity. At this moment, no new constraints are generated, and we have 
to compute the last Poisson bracket { f d?x€H, f dynH}pp. This Poisson bracket is too 
complicated to analyze, instead, we consider a simple case that is \/g€ = 6°(x — y) with n 
an arbitrary scalar function. In the process of computation, the following two formula are 
used, 


1 
ôRij = 59 (Vi V 559K: + ViVidgnj — ViV kôgij — ViVjôgki), (2.22) 
i l; 
Vin = 59 "(Vj 69n + Vrôgjı — Vidg5n)- (2.23) 


After a tedious calculation, we obtain the following result expanded with respect to the 
covariant derivatives of 7 in different orders. 


1... 5: , 
{H(x), / dynH}pp = 2G (EVE Vin + BYV Vin + Vin twn), (2.24) 
W 
with aJ* given by 
atsk = (Cum gI% 4 Chim gs = Cuan _ Cg) Kim, (2.25) 


where CJ is defined as enc E, in which C,* = gmC™ with C™ the Cotton tensor 
defined before. The term including the second order covariant derivative of 7 takes the 
form 


BIV;Vin = VGVinV y Kimo g* — KimC™*g™), (2.26) 


where the notation (kji). = (ijk + jki + kij) represents the cyclic permutation among 
indices i,j,k. The term proportional to the first order derivative of 7 is a little lengthy, 
and we have to introduce some abbreviations 


pidklm = Citim gk ra Cid ght pi Gim gik (2.27) 
gidklm _ Cum okt i Cid ght E CIM gik, (2.28) 
Then 
Vin = FIT GNV mV o) Kj + KyiV GNV kV mt 
+28" FEV nV iV pK jm + 2K jmVinV pV gs 
oO" Ri + C'i RI + C" RY) KaVin 


1. A a ii 
a5 jn) — Ri CO) KVin (2.29) 
The last one w is 


w = Vi(C™ RI Kj + COR Kay + OË RK) 
+CU" (VVVK + ViVi Vi Ky — ViV' Vik jk — ViVkV jK) 
HK! VVV: + KY, Vj ViVi — Kijk VVV: — KVVV) C. (2.30) 


At first sight, one may think that constraints from the third covariant derivatives of 7 
are aJk=(, it is not true. The reason is the following. We first notice that 


VeViVin = VkViV jn. (2.31) 


This originates from the formula of second order covariant derivatives of a scalar function, 
which is 

V;Vin = ðjðin — TE ôn, (2.32) 
where r$, is the Christopher symbol. For a torsionless space, r$, is symmetric in the 


two lower indices, combining with commutativity of partial derivatives, we have Eq. (2.31). 
Furthermore, we have the following identities, 


1 1 1 
VaVjVin = 5 V&V Vaen + (Vk ViVi — ViVe V5) + 5 (VeViVi — VaVe Vi)n, (2.33) 


where we used Eq. (-31). With the help of Eq.(2.31), we find that the first term in above 
equation is completely symmetric in three indices. If we denote the symmetrization among 
three indices by (ijk) = 4 (ijk+jki+kij+jik+kji+ikj), then V&V; Vien = V&V; Van. 
The remaining two terms in Eq. (2.33) can be reduced to the first order derivative of 7 by 
utilizing the following formula, 


(VE VaV = ViVi ViN = Ri aN: (2.34) 


After this process, it is clear that the third order covariant derivative term in Eq. (2.24) 
can be written as 


ET V5 Van: (2.35) 


Then inheriting the symmetry of V(,.VjVj)7, the effective components of aJ® are (3), 
Explicitly, they take the following form 


alik) = OR iy CIm gik $ CPP GI K mi 7 5 (Gig 4 Ori gim a CF g ™) K ma 
1 psy. Papa a As . 
= (Cree + Cragin + OIE G) K ma. (2.36) 


Because covariant derivatives of 7 of different order are independent, consistency requires 
that the coefficients in front of covariant derivatives of 7 should vanish (Because H is 
a first class constraint and the reason will be given later). We deduce from the term 
al FV Vj Vn that 

atik) — 0. (2.37) 


To see whether this condition gives rise to new constraints, we will work in a special frame 
where Eq.(2.37]) becomes simple. First we write the 3-dimensional space metric tensor 
gijdz’dz’ in terms of vielbein, namely 


gijdx'da = 540°0°, (2.38) 


where ôa is the Kronecker delta function, new basis 6% is related to dx’ by 6% = Olda’. 
Since Cotton tensor is a symmetric tensor, we can diagonalize it through an orthogonal 
transformation O (x) at each given point. So in terms of new basis defined by ĝe = On Gg 
the component of Cotton tensor becomes 


coe = (OCOT) = diag{C}, C2, = Cı — C2}, (2.39) 


where Cy, Co,and — C1 — Cho are three eigenvalues of Cotton tensor, where we used the 
property that Cotton tensor is traceless. Contracting two sides of Eq. (2.4), we obtain 


T au 
Jgk = + (m = g” Tij), (2.40) 
this helps us to reexpress Kj; by nti 
. w aoe 
VJK” =T” — zI T (2.41) 


Now, we find that Eq.(2.37) gives the following seven independent equations 


Cin? =0, Cor”? =0, (2.42) 
Cır” =0, Cor! = 0, (2.43) 
Cir? =0, Car” = 0, (2.44) 
(Ci + 2C2)r!t = (20) + Co)n”? + (C1 = Co)r? = 0. (2.45) 


The solutions of these constraints can be separated into two classes 


e 1) Ci = C2 = 0 or in other words Cotton tensor should vanish. We notice that the 
other coefficient tensors also vanish because they all consist of Cotton tensor and its 
covariant derivatives. Since the fact that a tensor vanishes is a coordinate independent 
statement, these constraints are still required to be satisfied, even after one uses the 
three coordinate transformations to fix three components of gij. In a special frame, 
we see that the vanishing of Cotton tensor gives two constraint equations, actually, 
in a general coordinate frame, the vanishing of Cotton tensor indeed provides only 
two constraints. The reason is that the following properties of Cotton tensor 


Cij =C, 9 Cig =0, VjC” =0, (2.46) 


make Cotton tensor itself have only two independent components. Gauge-fixing and 
vanishing of Cotton tensor altogether make only one degree of freedom in gij be 
physical. Now, the Hamiltonian constraint and three momentum constraints elimi- 
nate four conjugate momenta, leaving two components of 7’ be physical. Altogether, 
the phase space is described by three unpaired fields. 


° 2) 
C1 #0 or C2#0, n? =0, n =0, cB =0. (2.47) 


This case is rather bad, it indicates that all the conjugate momenta of gij are unphys- 
ical, since three momentum constraints already elminates three of the six conjugate 
momenta. 


To interpret above results with the approach given by Dirac , we consider the time derivative 


of the constraint. Utilizing Eqs.(2.20), (2.21}) and (2.24), we obtain 
dH; 


P7 = NH + (OINİYH; F 0;(N7H;), (2.48) 
Ct = (ONY + AN, (2.49) 


where the operator A is defined as 
1... 2 
A= —2V age (0 AVANGV + BEN 5Vi +7'V;i tw), (2.50) 


where the coefficients aJ*, B, yê and w take the same expression as in Eq.(2.24). The 


preservation of constraint in time require ari ~ 0, aH ~ 0 on the constrained phase 


space. The first one is satisfied because H ~ 0 and H; ~ 0. While to satisfy the second 
one gives a differential equation AN ~ 0. It is remarkable that A has no inverse on 
the whole constrained phase space, due to the existence of configuration gij, mI making 
the coefficients in front of covariant derivative vanish. As we have analyzed before, the 
configuration gi; with Cj; = 0 can achieve this. Therefore, H cannot be perceived as a 
second class constraint, since to define the Dirac bracket associated with the second class 
constraint, the inverse of A is indispensable. So H can only be a first class constraint, and 
its Poisson bracket generates constraints according to the property of first class constraint. 
This justifies our previous treatment of the terms yielded by {H(x), f d'ynH}pv. 


To summarize, we have found new constraints generated from the Poisson brackets 
of Hamiltonian constraint. These new constraints reduce further the phase space in a 
way that it appears that no symplectic structure exists, or eliminate all the degrees of 
freedom. We expect that the Poisson brackets among the new constraints and H, Hi yield 
more constraints, until the constraints form a closed algebra. When this is done, all the 
constraints are called the first class. Note that the Hamiltonian density is not a second 
class constraint, since the Poisson bracket obviously does not have an inverse. 

Put together, all new constraints either eliminate all degrees of freedom, or make the 
reduced phase space unphysical. 

A complete discussion deserves to be carried out in another work. 


3. Discussion and conclusion 


In previous section, our calculation shows explicitly that H and H; do not form a closed 
algebra as what happens in general relativity. Intuitively, we feel that non-closure of 
the Poisson brackets among H and H; can be interpreted by considering the relationship 
between constraints and gauged symmetry . The details are what follows. One can calculate 
the Poisson bracket between the combined constraint NH + NH; and gij then obtains 


g d?xn(NH +N'H;i), gijtpb = —OinN; — O;nNi — ngiz. (3.1) 


It is nothing but the variation of gi; under (3+1)-dimensional coordinate transformation 
in time direction. So the four generators of (3+1)-dimensional coordinate transformation 
are (NH + N'H;) and H;. This seems to suggest that the constraints H and H; require 
the theory to have full (3+1)-dimensional covariance to respect them, or more constraints 
should be added. The reason is that in a field theory, a local constraint is always accom- 
panied by a gauge symmetry. if we denote the constraint by C, the meaning of constraint 
can be expressed by 

C|phys) = 0. (3.2) 


This equation implies that the physical state is gauge invariant if there is a corresponding 
symmetry. In H-L theory, although there is no diffeomorphism invariance in the time 
direction, Eq. (B) still looks like such a transformation and imposing this constraint is 
somewhat in conflicts of the starting point?. 

What we have discussed previously is based on the UV action of H-L theory. To 
complete our discussion, the full action containing the description of H-L theory both in 
UV and IR region should be taken into account. In the IR region, some operators with 
lower mass dimension will become relevant. The detailed balance principle forces the action 
to take the following form Bi 


SS J dxan KK" -— AK?) — pao 


?Imposing local Hamiltonian constraint, but lacking the diffeomorphism invariance in the time direction 
may cause strong coupling problem A in the IR region of H-L Theory. Before the revised version of our 
paper appears, paper ip) also found this independently. 


2 oa 
+e Ra VR} " E RyR” 
W 
E u g — 4A 
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R? + AwR — 347 )}, (3.3) 


For this action to be a deformation of Einstein-Hilbert action in the IR region, it is natural 
to express this action in relativistic coordinates by rescaling t, 


£o = Ct, (3.4) 


with the emergent speed of light and effective cosmological constant given by 


_ |_^w = 
c= 3V TL3 A= JAW: (3.5) 


At the same time, parameter À should be equal to 1 for the kinetic term taking the same 
form as its counterpart in general relativity. Thus to have a real speed of light, Aw should 
be negative. Setting A = 1, the action becomes 


ij 1 P 
S= J uds aN (Ky K" =R- cue" 
2 
H ij H ij 
tage V - E Rg R” 
p3 
+o (GR -AwR + 3Aiy)}, (3.6) 
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The momentum constraints corresponding to this action remain the same form as in 
Eq.(2.19), the Hamiltonian density contains more terms than in Eq.(P.11)). We denote 
the new Hamiltonian constraint by H to distinguish it from the previous Hamiltonian 
constraint 


y l ;; 1 . 7 2 7 
HiV/9 = gn Gaer” —R+2A+ ga Cue” = aa €I Ra V; Rh ia E RyR” 
W W 


——R?, (3.7) 


where the speed of light has been set to 1 and then the relation Aw = —8/? has been 
used. It is noticed that the first three terms make up of Hamiltonian constraint in general 
relativity. Calculation of Poisson bracket of H is similar to the previous one, and we obtain 
the following result, 


2 ee i e n 7 
{H(x), | Pyn = -2vgzr (a IRD. Vin + B VVN +Y Vin + Gn), 
W 
+8;(ng Hi) — g Hið;n. (3.8) 


The last two terms come from the well known result of Poisson bracket between the Hamil- 
tonian constraint of general relativity. In the previous section, our discussion mainly con- 
centrates on coefficient tensor of the third order covariant derivatives of 7, we want to see 


whether this new coefficient tensor will change our results. Then we compute a and 
find that it amounts to replacing the Cotton tensor C” in attt by CY — Mw pis . The 


new contribution -EW RY comes from term = tE Ra V R}, since it also contains the 
WwW 


Pp 2 b 55 
third order derivative of metric gij. If we denote the tensor C% — Mw pij by &”, then the 


vanishing of @3*) becomes the following seven independent constraints, 


(£i — a)r? = 0, (£2 — tein = 0, (3.9) 
(£1 — X2) = 0, (Xp — X3)r" = 0, (3.10) 
(Sy =¥2)a* = 0, (£2 — X3)r” = 0, (3.11) 
(£2 — ela + (£z — ar + (£1 — a = 0. (3.12) 


To obtain above equations, we also work in basis where X;;j is diagonalized by 44; = 
diag{d,, £2, 3}. the three eigenvalues are independent because usually 5” does not 
satisfy the traceless condition. Despite this difference, our analysis can still be applied to 
this case and our conclusion is unchanged, the same difficulty remains. 

One may choose not to impose the Hamiltonian constraint, this contradicts the require- 
ment that the lapse function is a full-fledged function. Moreover, if one does not impose 
this constraint at the beginning, how can one obtain the usual Hamiltonian constraint in 
Einstein theory in the infrared regime? (This constraint is crucial in going back to the 
Newtonian limit) Finally, we want to remark that even if one can come up with some cure 
of this problem, it will be very difficult to come up with modified theory containing a spin 
2 graviton. 

Note added: When this paper was reviewed, we were informed that another paper 
holds a different opinion on the constraint structure in Horava-Lifshitz theory. We 
disagree with the result in (id) for the following reason. Most of discussion in that work 
is based on a perturbative method, while this method is not suitable for discussing the 
fundamental degrees of freedom of a theory. For instance, the new constraints obtained 
in our paper are at least of the second order perturbations around Minkowski background 
(adopted in {10]), so at the linear level they do not show up. However, these constraints 
already determine how many degrees of freedom a theory can have before one carries out a 
perturbative calculation. In summary, our point of view is that this perturbative method 
may be useful in solving equations of motion but is invalid for counting the number of 
degrees of freedom of a theory. 
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